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SUMMARY 


This report documents research performed under Contract NAS 8-21082 
between 1 January 1967 and 31 October 1967* The object of this program 
is to measure the absorption coefficients and fine structure parameters 
of water vapor at wavelengths between 1 and 10 (j. and at temperatures be- 
tween 1200° and 3000°K and to incorporate these into a practical analytic 
procedure which may be used for the evaluation of the radiative base heat- 
ing of Saturn V and other launch vehicles. Scientific and technical per- 
sonnel who spent appreciable amounts of their time on this study were 
C. B. Ludwig (Principal Investigator and Project Leader), J. A. L. Thomson, 

W. Ife-lkmus, M. L. Streiff, C. N. Abeyta, R. Janda, and D. Suttie. The 
Contract Monitor was R. M. Huffaker, Aerodynamics Branch, National Aeronautics 
and Space Administration, Jfershall Space Flight Center, Huntsville, Alabama. 

A uniform volume of high temperature water vapor was produced above 
the combustion zone of a long three inch wide flat burner which was con- 
structed for these measurements. Emission and absorption measurements 
were made with a pathlength of 20 feet. 

The results obtained under the present program have been combined 
with previous data in order to generate tables of water vapor absorption 
coefficients covering the spectral region 50 to 9300 cm - '*' and the tempera- 
ture interval 300° to 3000°K. The representation of the inverse of the 
line spacing covers the intervals 1150 to 7500 cm -1 in the temperature 
range from 600° to 3000 °K. 

Several nonisothermal measurements with four five-foot sections at 
different temperatures were made to test the nonisothermal radiance calcu- 
lation procedure. 
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A theoretical study of the sensitivity of band models to deviations 
of the line shape from the assumed Lo rent z shape was made and plausible 
line intensity distribution functions were, investigated. 
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SECTION 1 


j 

INTRODUCTION j 

Measurements of hot water vapor absorption up to a pathlength of ten 
feet were carried out previously. Assuming the applicability of the 
statistical band model,- values for the absorption coefficients and fine 
structure parameters were obtained. 

The purpose of the research reported here is the extension of these 
measurements to twenty feet, in order to obtain more reliable values of 
the fine structure parameter in spectral regions where the square -root 
region of water vapor becomes better defined only at longer pathlengths . 

These spectral regions include the troughs between the vibration-rotation 
bands and the higher order bands at elevated temperatures . 

Of course, at and near the bottom of the troughs, much longer path- 
lengths than the present twenty feet are needed to bring the water absorp- 
tion into the square-root regions. However, the present results are suf- 
ficient as to extrapolate with confidence the values of the fine structure 
parameter into these trough regions, thereby eliminating the uncertainty 
existing in the previously reported values . 

New values of the absorption coefficients k and the reciprocal of the 
average line spacing d are reported. The decision to present l/d rather 
than "a" was dictated by the desire to reduce the complexity in the calcu- 
lation of e at other than the present conditions . The representation of the 
line half width y as a function of self and foreign gas broadening, tempera- 
ture, and pressure is the same as given in Refs. 1 and 2. 

The new values of k are very similar to the values reported previously. 

The new values of the fine structure parameter a (and therefore l/d = a/y) 
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are similar to the previous band-averaged values of a for the 6.3- and 2.7-d 
band at all temperatures, for the 1.8- and l.l4-|j, band only at temperatures 
below 1500 °K. Higher values of a are obtained for these bands at higher 
temperatures (> 1500°K) and in the troughs. In order to facilitate the 
computational procedures in complex programs, it is desirable to represent 
the frequency dependent term l/d by simple equations. Previously, band- 
averaged values were suggested. We propose now to replace the band-averaged 
values with a sinusoidal function, which simulates much more closely the 
measured and extrapolated values of l/d. While this does not render the 
computations more difficult, it gives better values of e at long and non- 
homogeneous pathlengths than by using band-averaged values . 

In addition to these measurements, the independent control of four 

separate five-feet sections was utilized to make inhomogeneous measurements . 

The comparison with calculations using the present set of k and l/d and 

the Curtis-Godson approximation gives good agreement, even for the extreme 

case, where the radiation of 5 feet of H^O at 2500 °K is passing through 15 

feet of HgO at 1200 °K. Although we estimate that the calculation of the 

spectral and total emissivities of water vapor above 1200 °K is accurate to 

within ± 20 $ for homogeneous and inhomogeneous paths at most of the conditions 

relevant to the exhaust of the Saturn V, certain additional studies must be 

made to verify this confidence level to temperatures below 1200 °K and to 

more extreme inhomogeneous paths . The first task can be carried out through 

3 

the analysis of the extensive experimental data taken by Rocketdyne and by 

4 

the University of Michigan. The second task requires a combination of 
experimental and theoretical investigations. Mere extreme inhomogeneous 
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path measurements must be taken and work on different models for use in 
the Curtis-Godson approximation must be continued. 
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SECTION 2 

EXPERIMENTAL PROCEDURE 

The installation of the ten-foot burner was described previously. 1 For 
the present program, ten feet of burner with the necessary propellant, 
cooling and control system were added. A photograph of the completed 
burner inside the tank is shown in Fig. 1. At both sides of the burner, 
the manifolds connecting five one-foot sections can be seen. The four 
infrared safety devices (one for each five -foot section) are located 
outside the tank and cannot be seen in the figure, except one hole in 
the tank wall near the entrance, through which the collimated radiation 
passes. The infrared safety devices consist of a FbS detector, an 
amplifier, and a relay, which governs the two main control valves for 
the propellants . If the radiation drops below a previous set level, the 
main control valves are shut and a nitrogen purge is initiated. In ad- 
dition, the output of these amplifiers is displayed on a Sanborn recorder 
in order to check the constancy of the radiation during a spectral scan. 

The maximum variation of the output level was within the accuracy of the 
control gages (± 2$>) . The check of the homogeneity was again made with 
a thermocouple, traveling the total length along the line of sight. The 
RMS value of the variations was less than 5$. 

A total of 2b runs were made. In order to establish consistency with 
the previous runs, six runs with a five-foot pathlength using different five- 
foot sections, six runs with a ten-foot pathlength, and twelve runs with a 
twenty-foot pathlength were made. During the course of these runs, one 
error was discovered, which existed in some of the previous runs made last 
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year at ten-foot pathlengths. A discussion of this is given in the next 
section. 
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SECTION 3 


DATA ANALYSIS 

The same data reduction procedure as described in Ref. 1 was used. 

Previous runs at 2' , 5' , and 10' together with new runs at 5 ' , 10 ' , and 

-1 

20' were used. The spectral emissivities were averaged over 25 cm 
intervals . The averaged emissivities were then least-square fitted to 
2nd degree polynomials as functions of the temperature. The RMS of the 
deviation to these fits was about the some as was found previously, less 
than ± 5# in the bands, but somewhat greater in the troughs, where the 
signal-to -noise ratio was smaller. However, the errors in the 20' trough 
measurements were less than the corresponding measurements at smaller 
pathlengths . 

From the curve-fitted data, the emissivities were calculated at 
specified temperatures between 1250° and 2750 °K. Using the statistical 
model with exponential line intensity distribution, the parameters k and 
a were obtained from the curves of growth at each cu. The curves of growth 
£n(l/l-€ ) were transformed into linear functions of u: 

^u/0n (l/l-e)) 2 = k -2 + u(Uak)” 1 . (l) 

From the intercept with the ordinate and from the slope, the values of k 
and a, respectively, are obtained. 

In comparing these data with the previous ones, given in Ref. 1, the 
following observations are made. 
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1 . 


The absorption coefficients are essentially the same, except 
in the troughs, 

2. The fine structure parameters "a" are larger for the higher 
order bands and for the trough regions between the bands . 

The reason for the difference is twofold: first, the 20' 

burner data improved the determination of the square -root 
region, which was previously well-defined only for the regions 
having the largest values of k; second, some of the previous 
10' burner data were in error by about 20$. 

This error occurred through a wrong calibration in the propellant 
system which influenced the determination of the temperature. After the 
data were corrected and new 10' data obtained as a check, they were then 
used in the present data reduction, which subsequently gave higher values 
of a. As an example, the curve of growth at 2000 °K for uu = 3500 cm ^ is 
plotted in Fig. 2. The previous data are given by (a), and the resulting 
curve of growth (least-square fit) is given by the dashed line. The present 
data are given by (A ) and the curve of growth is given by the solid line . 

The present curve has a smaller slope and, therefore, a greater value of 
"a" for the same k. At the same time, the a (RMS of the scatter of data 
points) is much reduced. 

As in the previous data reduction, where the largest pathlength was 10', 
the present data up to 20' are still not sufficient to determine the fine 
structure parameter in the troughs between the bands . These regions were 
presented in Ref. 1 as dotted lines because of this uncertainty. However, 
a certain consistency in the data let us believe that "a" was small in these 
regions. We are now convinced that the consistency was fortuitous. The 
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present data indicate a rise of "a" toward the wings of the bands so that 
the actual minima of "a" occur in the center of the bands (excluding the 
Q -branch) and not in the troughs. 

Although the absolute levels of "a” and the dependence of "a" on the 
wavelength are different as compared with the results given in Ref. 1, the 
difference in the total emissivity based on the two different sets of "a" is 
small and is within the uncertainty limits of the band model parameters. In 
Rigs. 3 and 3a, the total emissivities between 1150 and 5900 cm ^ at 2500° 
and 1500 °K, respectively, are computed based on previous k and band-averaged 
a's and on present k and l/d at a total pressure of 0.1 atm and 50# H^O and 
50# Wg mixture. The abscissa is the pathlength from 0.1 to Kr cm-atm. The 
maximum deviation of about 25# occurs at around 100 cm-atm for 1500°K. 
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SECTION 14- 


RESULTS 

In this section, we present the results of the absorption coefficients 
and the fine structure parameter in graphical and tabular form. Comparisons 
between experimental and calculated data, error considerations and results 
of inhomogeneous data are also given. 

4.1 THE ABSORPTION COEFFICIENTS OF H 2 0 

Data were obtained in the temperature range from 1250° to 2750 °K. Extra- 
polations were made to 1000° and 3000 °K, Previous data at 300° and 600°K 
were re-examined and changes incorporated where necessary. Comparisons of 
the new data versus previous data'*' are shown in Figs . 4-10 at the tempera- 
tures of 300°, 600°, 1000°, 1500°, 2000°, 2500°, and 3000°K. The major 
differences occur in the troughs only. The change of the trough values at 
4600 cm”'*' for 1000°K comes through an error in the previous data. No experi- 
mental data exist for the trough values at the temperatures of 600° and 300 °K. 
The new data are logical extrapolations from the high temperature results . The 
results are given also in tabular form in Table I. 

4.2 THE INVERSE OF THE LINE SPACINGS 

While the format and the actual values of the absorption coefficients are 
changed very little, the representation of the fine structure parameter has 
been changed. In the case for pure collision broadening for E^0 spectral lines, 
we write 
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where 


■£° ■ ZK’^„ *, (9 *'•] * (*•$„ v (9 ■ o) 

J 


Preliminary data of the individual line half widths for different broadeners 
are the following (the same as given in Ref. 2). The data in parentheses 
are guesses. 

* H 2° 

y h 2 o= ° M 


h 2 0 

^0 


= ( 0 . 


y 


= 0.09 


h 2 0 


= 0.04 


HpO 


= (0.05) 


y 

r co„ 


= 0.12 


y 

r co 


= ( 0 . 10 ) 


Values of the inverse of the line s pacings {i/d') are given in this report. 
These are identical to the (l/d Tt) ), used in Model 3a. For Model 3. values 

J_it\ 

of l/d- o are used, where 


_1_ l_ 1+e 

d T „ d .. 

LR o 1-e 


-6/2T 


•072T 


( 4 ) 
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with 9 = 230CTK. For a binary mixture, Eq. (2) reduces to 


£° - P T [<0.44 (*E) ♦ 0.09 (^) 1/2 ) ♦ (i-=) (^i) 1/2 v“ 2 °] 


1/2 H 
*3 


(5) 


HO HO 

which has been used (with y . = Y n = 0.04) to reduce our experimental 

J t/g 

data of "a" to the values of l/d (- l/d^). The concentration c for our 

conditions was calculated from a the rmochemi c al equilibrium model. 

The reduction of the experimental values of "a" to the parameter l/d 

assumes the validity of the representation of the collision half width 

* H 2° 

through Eq. (3) or Eq. (5) and the three individual half widths y„ 0 , 

H 2 0 H 2 0 H 2 0 tt 2 U 

Y h q, and Yq • If In the future a new representation of Y c and/or 

new individual half widths are found to be more applicable, the fine 

structure parameter "a" can be obtained through Eq. (5) and the values 

of l/d listed in this report. The concentration e is given by 


c = -.1002 + .2802 X 10 -3 T 1 -.1089 X 10" 6 T 2 + .0291 X 10 -9 T 3 (6) 


Data were obtained for the temperature range from 1250° to 2750 °K in the 
approximate spectral regions : 


II50 - 2200 cm 
3200 - 4100 cm 
4750 - 5700 cm 
6700 - 7500 cm' 
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These spectral regions correspond to the major portions of the four 
vibration-rotation bands of H^O. In all cases, the values of l/d increase 
toward the wings of the bands and become undetermined in the center of the 
troughs and the values of l/d were extrapolated into the trough regions. 
Since the plots of the experimentally determined values of l/d versus 
temperatures were smooth curves (see Figs. 11-20), they were extrapolated 
to 600°K and to 3000°K. A verification of these extrapolations to 600°K 
will be made in the future through the use of Simmons and Rocketdyne data. 

The set of l/d versus u) for T = 600°, 1000°, 1500°, 2000°, 2500°, and 
3000°K is shown in Fig. 21, where all regions, which were extrapolated by 
us, are given as dashed lines. A list of the values is given in Table II. 

In order to reduce the number of values to be stored in a computer 
program, the l/d were represented by a sinusoidal function: 

Bn l/d = A sin (Bcd+C) + D(T) (7) 

where 

A = .7941 

B = .0036 
C = -8.043 

D(T) = -2.295 + .3004 x 10 -2 T -.366 X lo " 6 T 2 
In order to check this approximation, integrated emissivities were calcu- 
lated based on (l/d) from Table II and on (l/d) from Eq. (7). The calcu- 
lations were performed for a gas mixture of 50$ HgO and 50$ at a total 
pressure of 0.1 atm. The optical path ranged from 0.1 cm-atm to 10 ^ cm-atm. 
In Fig. 22, the total emissivities integrated from 1150 to 7500 cm ^ are 
compared. In Figs. 23 and 24, the band emissivities at 1500° and 2500°K, 
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respectively, are compared. The difference between the results is small. 

4.3 COMPARISON WITH EXPERIMENTAL DATA 

The consistency of our results was checked by comparing the experimentally 
measured emissivities with calculated values using the k and l/d from Tables 
I and II. The results are shown to be consistent (see Figs. 25-32) for all 
pathlengths (2, 5> 10, and 20 feet) at two representative temperatures 
(1500° and 2500 °K). 

As an independent check, comparisons with Simmons’ data are made for 
temperatures greater than 1000 °K. The results are given in Figs. 33 and 34 . 

4.4 ERROR LIMITS 

It is difficult to assign specific error limits to the data we presented, 
because they depend on many different parameters . We can estimate the 
precision with which the data were taken. From the knowledge of the 
statistical scatter, the accuracy of the intermediate and the final results 
can be estimated. It was found* that both precision and accuracy depend 
upon the spectral regions, the temperature and the pathlengths. In addition 
to these considerations, the uncertainty of the models themselves must be 
assessed. 

In order to give some overall error limits, we estimate the calculation 
of integrated radiancies (or total engineering emissivities) based on the 
data presented here to be within ± 5 % for pL 5 20 ft-atm. For greater path- 
lengths, the uncertainty is surely greater because it depends then more on 
the individual foreign gas broadener parameters, which are uncertain and, in 
some cases, are only guesses. For the calculation of spectral radiancies 
(or spectral emissivities), we estimate the overall error limits to be within 
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± 20$. There are several exceptions to this value. The center of the 
2.7-h band at the higher temperatures is known better than ± 10$, while 
all the trough regions and bands at the lower temperatures are probably 
known with less certainty than ± 20$. 

We will now discuss the individual contributions to the overall 
error limits. The precision of the individual data (in p.V) is largely 
determined by the signal- to -noise ratios of the system, the knowledge of 
the calibration blackbody source (energy, temperature) and the spectral 
calibration (micrometer drum divisions versus wavelength). The signal- 
to-noise ratio changes from 1 for low radiation levels to 100 for high 
radiation levels . Low radiation levels exist in the bottom of the troughs 
and in spectral region > 7000 cm”^ at temperatures < 1500°K, and no reliable 
data are obtained in these regions. The precision in determining the black- 
body temperature is within ± 1°K. The spectral calibration is precise within 
± 3 cm”^. 

The accuracy of the spectral data depends upon the conversion of the 
output of the amplifier (p-V) to emissivities . In the case of absorption 
spectra, the accuracy is determined by the precision with which the data 
were taken, and is, in general, poorer by a factor of 2. In the case of 
emission spectra, the precision with which the gas temperature is known, 
in addition to the precision of the amplifier output, determines the accuracy 
of the spectral emissivities. The temperature of the gas is known within 
± 3$, This value was established through the precision of the flow meters 
and through the statistical scatter of many temperature determinations . The 
± 3$ precision in temperature means an inaccuracy of ± 5$ to ± 30$ in the 
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spectral 6missivitl.es from 3000 cm ^ at 3000 °K to 8000 cm ^ at 1000 °K. 

(See Pig. 35). 

An improvement in precision is effected by the curve fitting process 
for e versus T. We found that the RMS value of the deviations of the experi- 
mental spectral emissivities is ^ ± 5$ for the major portions of the bands 
and ^ ± 20 °/o in the troughs down to the lowest temperature of 1250 °K. 

The deduction of the absorption coefficients and fine structure 
parameters from the experimental data depends upon the knowledge of the 
state parameters of the gas, the pathlength, and the band model employed. 

As stated before, the temperature is known to within ± yjo; the total 
pressure was always assumed to be 1 atm; the partial pressure for water 
vapor and oxygen was calculated based on the assumption of chemical 
equilibrium, and the uncertainty is unknown; the geometric length is 
uncertain due to the temperature gradient at the two ends of the burner 
in the zones where the hot gas mixes with the cold nitrogen flow. Thermo- 
couple readings at low temperatures and visual inspection at high tempera- 
tures, where the flame exhibits an orange color, indicates a transition 
zone thickness of approximately 4" . At a nominal pathlength of 20 feet, 
these 4" amount to an error of about 2$. A value for the uncertainty in 
using a particular band model has not been established as yet, but will 
be done in the near future. A measure of the validity of the statistical 
model can be seen in the fit of ^u /S/n. versus u to a straight line 

[see Eq. (l)] . Figures 36-38 show the RMS scatter (denoted by a) in the 
fit of the experimental data to Eq. (l). This is not readily interpreted 
in terms of accuracy of k and a. If the measurements are primarily near 
the linear region of the curve of growth, the points in the transformed 



equation [Eq. (l)] have nearly equal ordinates. Thus k is well-determined, 
with experimental accuracy of about § a. The slope of Eq. (l) may be poorly 
determined, and the error in a may be large. If the measurements are primar- 
ily near the square-root region, the slope of Eq. (2) may be well-determined, 
but not the intercept. Thus, the product ka is known, but k (and, hence, a) 
may not be known separately with, any accuracy. If the measurements are 
located primarily along the transitional region of the curve of growth, k 
and a may both be known with comparable accuracy, but with error greater 
than ■§• ct. Only if the measurements extend well into both the linear and 
square-root regions (so that c represents the uncertainty in k as well as 
in the slope ka) will both k and a be determined with an accuracy of about 
\ cr and o, respectively. The trend exhibited in Figs. 36-38 is to be 
expected. Low o values exist in the bands and at the higher temperatures, 
while higher values exist in the troughs and at the low temperature. 


4.5 INHOMDGENEOUS GAS MEASUREMENTS 

In addition to the homogeneous measurements, the independent control of 
four separate five-foot sections was utilized to make inhomogeneous measure- 
ments . The experimental results are given in Figs . 39-44 in terms of 
/ 2 <-l 

W/cm cm for the following conditions. 


Section 

A 

B 

C 

D 


Temp 


Cone 

Fig. 39 

Fig. 40 

Fig. 4l 

T 

2480 

1210 

1210 

c 

•37 

.1.4 

.14 

T 

1220 

1520 

1210 

c 

.14 

.18 

.14 

T 

ll8o 

2000 

ll8o 

c 

.13 

.26 

.13 

T 

2480 

2480 

2480 

c 

.37 

.37 

•37 
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The experimental data are compared with calculated ones, using the Curtis - 
Godson approximation with Madel 3 and 3a. The agreement of the theoretical 
and experimental values is within the uncertainty of the temperatures . The 
comparisons are made at 1 atm and the difference of the theoretical values 
determined by the two models is at most 12 $. In addition to the calculations 
at 1 atm, calculations were performed at 0.1 atm, where no experimental data 
exist. In that case, the largest difference of the theoretical values 
amounted to 50 $ between the two models, Jfodel 3a always being higher and, 
therefore, being conservative. 

In conclusion, we remark that the simpler model 3a is preferable, giving 
good agreement with experimental data at 1 atm and conservative values at 


lower pressures 



SECTION 5 

ANALYTICAL STUDIES 

Studies have continued on optimizing the selection of band models for 
use in particular situations . Certain results are summarized here in 
regard to two topics: (a) the effect on non-Lorentz line shape and 

(b) the choice of line intensity distribution function. 

5.1 NON-LORENTZ LINE SHAPE 

Experimental studies have been successful in demonstrating deviations 

from Lorentz line shape in certain situations, e.g., in the troughs 

5 

between the widely spaced lines of hydrogen halides, and in the region 

6 

beyond a sharp band head ( CO^ ) . In the former case, deviations were 

noted between 2 and 20 cm - "*" from the line centers, and in the latter up 

to 200 cm "*■ from the nearest line centers . 

In the case of HF, above -Lorentzian shapes were noted, and the wings 

were empirically fitted 5 to a form q'(u)-u)q)" , where n (which is 2 for 

lorentz lines) was observed to be of the order of 1.8. 

In the case of CO^ and CO, the line wings about 100 cm - "*" from the 

£ 

line centers were found to be below -Lorentzian in shape. In these cases, 
the experimenters assumed an empirical factor applied to the line shape, 
which factor decreased experimentally with distance from the line center. 

In studying the behavior of models of lines where the local properties 
do not vary abruptly, as near a band head, or in the far wings of a molec- 
ular band system where the local intensity may drop off exponentially, we 
are most concerned with the deviations in the near wings of the lines, 
since this may have the greatest effect on the local emissivity properties. 
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We will investigate the situation where the absorption coefficient in 


the wings of a line varies as (cjj-cjUq) - : 


k(u)) 


-n 


Ico-oDqI for |u)-(u 0 | large 


a) 


In analogy to the case of the Lorentz line shape, we assume 


k(u>) = 


, n | in ’ 

b + | (i)— (i) Q | 


( 2 ) 


where b has the usual meaning of half -width at half -intensity, and a is 
proportional to the line intensity. Prom the definition 


- J k(») 


duu 


( 3 ) 


we obtain directly by integration 


k (u,) = (Sb n ~Vn) [| n sin(TT/n)] 

v , n | in 

b + |u)-cu 0 | 

For n = 2, Eq . ( 4 ) reduces to the usual Lorentz equation : 


00 


k(cu) = 


Jb/lL 


v 2 ( \2 

b + (u)-oo o ) 


( 5 ) 


From the definition of equivalent width of an isolated line, 

00 

= J ^l-exp[-k(a)) u]j- dcu , 


W 


( 6 ) 


we have 


W 


- J —~f J} 


0 


b n + |uMU 0 l 


da) 


( 7 ) 


This equation is not immediately integrable . However, for investigating 
the strong line region, we may study the asymptotic form 
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w 


£ J { 1 - ex P C- au I W -<^ 0 l“ n ] | Clw, 


( 8 ) 


which yields 


1 1 
W * 2a 11 r (1 - J ) u n , 


(9) 


or explicitly. 


W * 2( | sin | ) n T (1 - i ) (ff ) Q (Su/h f h. 


n 


( 10 ) 


An approximate form of Eq. (10) can be obtained by expanding in 
powers of ( n-2 ) : 


| -(n-2)/4 

W S 2(Sbu) (Su/b) [ 1 + .05 (n-2) ]. (ll) 


For l.T ^ n£ 2.3, Eq.. (ll) approximates Eq.. ( 10 ) within 1$, so that Eq. 
(ll) still provides a very close upper limit. For n = 2, both Eq. (10) 
and Eq. (ll) yield 

1 

W £ 2(Sbu) 2 , (12) 

the usual equation for the "square -root” region of a Lorentz line in the 
strong-line approximation, i.e., when the line is strongly absorbed in the 
region near the line center. 

We note that in the strong- line region for a line of shape given by 
Eq. ( 2 ), we have 

1 

W °= u n . (13) 

• 455 

Thus for n — 1.8, we have W u in the strong line region. 
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If we consider a band composed of randomly located lines whose 
intensities are distributed according to some normalized intensity dis- 
tribution P(S), we have 


W = J W(S) P(S) dS 


(14) 


For the case of an exponential intensity distribution.' 

P(S) = S' 1 exp (-S/S Q ) , (15) 

and the k(co) given by Eq. ( 2 ), we have 


00 00 

f r 

f i c .™ f (^"Vff) (n/2) sin 

(ff/n)Su "I 1 

J S = 0 ‘ <0 ■= 0 

■s x — exp — 

> b n + I CO - CO 

|n J J 


xcko exp(-S/S^) dS . ( 16 ) 


After lengthy manipulations, including interchange of order of integration, 
we obtain an exact expression: 


1 


„ n 

W = S Q u [ 1 + n sin(ir/n) S Q u/ 2ffb] 


1 


(IT) 


For n - 


2, we obtain the customary expression 



[ 1 + S Q u/ ffb] 


1 

2 


( 18 ) 


On division by d, the mean line spacing^ Equation ( 17 ) yields, in 
the limit of small S^u/d, 


W w SqU 
d d 


J 


5-h 


(19) 



and, in the limit of large S Q u/d, 



If we identify these asymptotic regions with those of a model'' comprised 
of lines of an equivalent intensity Sg and equivalent spacing dg (using 
Eq. (lO)) we find 


n 




s 0 = S E [ ? 8in ( | 

) r (1 * 

n ) ) ^ n-1 

(21) 



a = a E[l >T( I- 

b]'^ 1 ■ 

(22) 

In terms 

of these parameters, we have 



W _ 
d 

V 

{l + [2sin(H)- 

n 

1 ‘ n - 1 





x V I 
2b J 

i-i 

n 

* 

(23) 


On introduction of the dimensionless parameters 

b e - arb/dj, 

and 

Xg = S E u/2ffb , 

Eq. (23) becomes 

d ~ ^E e n^ X E^ * 


( 2k ) 


(25) 


( 26 ) 
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where 


' n 
" n-1 


n 

n-1 


e n ( V “ *E { 1 + [ i? Bln < l > ] [ r f 1 - ; > ] 


i-1 

n 


Xirxg} 

Again, we note that for n = 2, Eq,. (27) becomes 

_ 1 

e 2 (x E> = *E [ 1 * I *E ] 


( 27 ) 


(28) 


the curve of growth previously designated"^ e(xg). 

In Fig. 42, the curves of growth e Q ( x g ) are plotted for n = 1.8 
(dashed), 2.0 (solid), and 2.2 (dotted). A wide enough range of is 

used to show the fusion of each curve with its asympt6tes. 
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5 .2 LIME INTENSITY DISTRIBUTION FUNCTIONS 

The selection of an appropriate line intensity distribution function 
should be based at least to some extent on the theoretically anticipated 
distribution of intensities in the spectral region which the model is 
to represent, since it would be extremely arduous, if not completely 
impossible, to make an actual count of lines over many orders of magnitude. 

Some models which are still commonly used presuppose intensity dis- 
tributions which are based at least as much upon mathematical expediency 
as upon physical reasoning. For instance, an intensity distribution in 
which all lines are assigned equal intensities is clearly not an optimum 
description for a high-temperature molecular spectrum. However, under 
this assumption, the curve of growth of a random band model is identical 
to that of a single isolated line. Similarly the assumption of an expon- 
* ential intensity distribution, while more reasonable, is not based on direct 
physical reasoning, but is employed because it removes the Bessel functions 
from the Iadenburg-Reiche function and leaves a simple algebraic form. 

The distribution of intensities within a hypothetical band is frequently 
described by a probability distribution P(s) which is normalized: 

DO 

P(S) dS = 1 (29) 

o 

Thus the probability that any arbitrary line has an intensity between 
and is given by 



The probability distribution function is proportional to the number density 
of lines with respect to intensity 

P(S) = (31) 
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Using the chain rule, we can write 


dn _ dn <3E 

dS " dE dS • 


( 32 ) 


The dominant factor in the general expression for line intensity is 
the Boltzmann factor. From the approximate relationship 


S « exp(-E/kT) , 


(33) 


we find 


dE q-'l 

dS * S 


(3^) 


As a result of effects such as nonrigidity and multiple vibrational degree 
of freedom, the number density of energy levels will -increase with E. If 
we postulate an approximate power law 


— « E m <x 

dE . 


m 

[V tys)] 


(ms o) 


(35) 


we have 


m 


P(S) = K S -1 [in (S M /S)] 


From the normalization requirement. 


P M 

I , p(s) 
V R 


dS = 1 


where S M is the maximum and S m /R the minimum intensity, we have 


M' 


K = (m+l)/(lnR) 


m+1 


(36) 


(37) 


(38) 
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We then have, explicitly. 


P(S) 


(m+1) S _1 [ln(S M /S)] / (lnE) 


m+1 


( 39 ) 


For m - 0, this expression reduces to 


P(S) = S -1 /lnR, 

which case has been discussed previously.^ 
The average value of S is given by 


S = J SP(S) ds = (m+l)(lnR)“ m " 1 J M £lnS M /s] dS 

V R 

which, for large R, is 


S = S m T (ra+2 )/(lnR) 


m+l 


Similarly, the average value of S 2 is given by 


(40) 


(41) 


(42) 


3 ” = Js 2 P(S) dS = (m+l)(lnR)~ m " 1 J M [lnS M /s] S~ 2 dS , (43) 

S m /R 


which, for large R, is 

3 ~” = ^ r (# 5 ) 2 , tl /( i ««)'* 1 • 


By defining the parameters and such that 


VS - s / a 


and 


( 45 ) 


«* 4 /«b = ^/ 4 


( 46 ) 
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we obtain 


a, - V^ 1 


( 47 ) 


and 


a*. - 


(48) 


The probability distribution functions for R = 10^ and m = 0,1,2, 
and 3 are shown in Fig. 43. Note that in the region of smallest intensities 
the slope is somewhat steeper than -1, If we define 


x = In (S M /S) 


( 49 ) 


and 


y ^ in P(S ) , 


(50) 


Eq. (39) becomes 


y - In £(m+l)/S^ (lnR) m+ ^ "j - x + m lnx , 


(51) 


from which we find 


- -1 + mx" 1 . (52) 

Thus for x oh -14 (S = S^/r) and m = 3> the slope « -1.2. 

It does not appear possible to obtain closed form expressions for the 
curves of growth associated with P(S) of the form of Eq. (39)* We can, 
however, note a P(S) which is proportional to S"^ (shown in Fig. 43 as a 
dashed curve, with slope -2) is far more heavily biased toward the lower 
intensities than the curves shown for m = 1,2,3. The normalized distribution 
is expressed by 
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(otherwise) 


(53) 


P(S) = S M (R-l )" 1 S ' 2 
= 0 

If this expression is convoluted with an exponential distribution 
(for convenience), which broadens the distribution of intensities even 
further the resulting distribution is given by 

PCS) = (R-l )" 1 S “ 2 [(S M + S) exp(-S/S M ) - (S M + RS) exp(-RS/S M )] . ( 54 ) 


By following the same procedure as in Ref. 7 , we obtain the following 
explicit expression for the curve of growth resulting from the P(s) of 
Eq. (54): 


W /i = 7 2 ( V< 


(55) 


where 




_ ^ 


InR 


In < 


2 4 

[l+2lTRX E /(lnR) 2 ]-l [l + 2ffXg/ ( InR) 2 ]* 


+1 


[ [l+^RXg/tlnR) 2 ] +1 [l+2JTXg/(lnR) 2 ]‘ 


(56) 


v-1 


and 0 g and Xg have the meanings defined previously. 

This curve of growth is shown in Fig. (44) as a dashed line. For 
comparison, the Ladenburg-Reiche curve (solid curve) corresponding to 
P(S) = 6 (S -Sg) and the curve ^ h(Xg) corresponding to P(s) = (S InR)" 
exp(-S/ffS_) (dotted curve) are also shown. The mutual asymptotes are shown 
as solid lines. 

The curve F^ (Xg) provides only a gross lower limit to the curves of 
growth which would correspond to the P(S) of Eq. (39).* shown in Fig. (43). 
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Fig. 8. Comparison of previous k (*) and present k (solid line) at 2000°K. 
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Fig* 21* Plot of (l/djj^) versus to between 1150 cm ^ and 7500 cm for T = 600% 1000% 1500 
2500 and 3000 °K. The extrapolated regions are indicated by dashed line s. 
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Fig. 22. Integrated emissivity between 1150 and 7500 cm 1 versus pL, based on (l/d) from Table II 
(solid line) and (l/d) as a sinusoidal function of a) (dashed line) for 1500° and 2500°K. 
Total pressure is 0.1 atm, 50$ BgO and 50$ N 2 . Note: The values for the 2500 °K case 

must be multiplied by 10. 
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of (jo (dashed line). Total pressure is 0.1 atm, 50 $ HgO and 50$ 
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Fig. 23 
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Fig. 2k. Integrated band emissivities at 2500°K versus pL, based on (l/d) 
from Table II (solid line) and on (l/d) as a sinusoidal function 
of <ju (dashed line). Total pressure is 0.1 atm, 50$ HgO and 50$ 
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Fig. 27. Comparison between experimental (solid line) and computed spectra (*) at 1500°K and 5 * pathlength. 
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Fig. 30. Comparison between experimental (solid line) and computed spectra (*) at 2500 °K a.nri 10 • pathlength. 
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Fig. 32. Comparison between experimental (solid line) and computed spectra (*) at 2500 °K and 20' pathlength. 
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Fig. 35. Error in e due to uncertainty of ± 3* in T versus u) for three different temperatures 
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Fig. 36 . RMS scatter of experimental values to curve of growth at T = 2750 °K. 
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Fig- 37- RMS scatter of experimental values to curve of growth at T a 2000 °K. 
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Pig. 38. ME scatter 
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Pig. 39 . Radiance from four slabs (X — 150 cm each) in watts/ cm 2 cm 1 . Solid line is the experimental curve 
Calculated values from Model 3 are given by (a), from Model 3a are given by (X). Conditions for 
T and c are given in Section 4.5. 




41/It/M SAN 



Radiance from four slabs {l « 150 cm each) in watts/cm 3 cm" # Solid line is the experimental curve 
Calculated values from Model 3 are given by (□), from Model 3a are given by (x). Conditions for 
T and c are given in Section 4.5. 
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Fig. 4l. Radiance from four slabs (jJ = 150 cm each) in watts/cnPcm” 1 . Solid line is the experimental curve 
Calculated values from Model 3 are given by (□), from Model 3a are given by (X). Conditions for 
T and c are given in Section 4.5 . 
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Fig. 43. Line intensity probability distribution functions. Solid line: 
P(S) « S" 1 . Dotted lines: P(S) * [0n (Sj/S)] m S _1 , for m = 1, 

2, 3» Dashed line: P(S) « S~ . All P(S) = 0 outside of region 
(Sj/R * s £ S M ). 
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ABSORPTION COEFFICIENTS OF H20 
K (PER CM AT STP) 
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AflSORPTION COEFFICIENTS of H20 
K (PER CP AT STP) 
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ABSORPTION COEFFICIENTS OF H20 
K (PER CM AT STP) 
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ABSORPTION COEFFICIENTS OF H20 
K (PER CM AT STP) 
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ABSORPTION COEFFICIENTS OF H20 
K ( PER CM AT STP ) 
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ABSORPTION COEFFICIENTS of H20 
K (PER CM AT STP) 

1/CM 300 K 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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ABSORPTION COEFFICIENTS OF H20 
K (PER CM AT STP) 

1/CM 300 K 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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ABSORPTION COEFFICIENTS OF H20 
K (PER CM AT STP) 

1/CM 300 K 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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ABSORPTION COEFFICIENTS of H20 
K (PER CM AT STP) 

1/CM 300 K 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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ABSORPTION COEFFICIENTS of H20 
K (PER CM AT STP) 
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ABSORPTION COEFFICIENTS OF H20 
K (PER CM AT STP) 

1/CM 300 K 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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FINE STRUCTURE PARAMETER (1/DLR) FOR H20 
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FINE STRUCTURE PARAMETER (1/DLR) FOR H20 
I/O 600 K 1000 K 1500 K 2000 K 2500 K 3000 K 
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FINE STRUCTURE PARAMETER (1/DLR) FOR H20 
1/CM 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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FINE STRUCTURE PARAMETER (1/DLR) FOR H20 
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FINE STRUCTURE PARAMETER (1/DLR) FOR H20 
1/CM 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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FINE STRUCTURE PARAMETER (1/OLR) FOR H20 
1/CM 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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FINE STRUCTURE PARAMETER (1/DLR) FOR H20 
1/CM 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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FINE STRUCTURE PARAMETER (1/QLR) FOR H20 
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FINE STRUCTURE PARAMETER ( I/DLR ) FOR H20 
1/CM 600 K 1000 K 1500 K 2000 K 2500 K 3000 
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